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As is well known, the existed perturbation theory can be applied to calculations of energy, state 
and transition probability in many quantum systems. However, there are diflterent paths and meth- 
ods to improve its calculation precision and efficiency in our view. According to an improved scheme 
of perturbation theory proposed by [An Min Wang, quant-ph/0602055 v7], we reconsider the tran- 
sition probability and perturbed energy for a Hydrogen atom in a constant magnetic field. We find 
the results obtained by using Wang's scheme are indeed more satisfying in the calculation precision 
and efficiency. Therefore, Wang's scheme can be thought of as a powerful tool in the perturbation 
\^ , calculation of quantum systems. 
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I. INTRODUCTION 



. The traditional method of perturbation theory tells us to calculate the energy and expansion coefficient step by step. 
That is, we should first calculate the zeroth order energy and wave function, and then the first,the second, and so on. 
In fact, it is easy to see that such a way introduces the approximation too early, and each order calculation of energy 
and expansion coefficient of wave function are based on the result of the former orders. After careful examination, we 
find that the traditional way does not consider the astringency of the expansion of the wave function. There is every 
possibility that if we reconsider the astringency, the result might be different. In fact, this is effectively embodied in 
y—^ ' Wang's scheme of perturbation theory [l]] . In this reference, the author did not introduce approximation until very 
T-H . late, and consider subtly and systemically the affection of high-order approximation to the low-order one. This finally 
^^O ' results in a different formalism of expression of perturbed solution of dynamics, and its expansion coefficients contain 
, reasonably the high-order energy amendment. 

■ It can been seen that some physics expressions were modified according to Ref. l] and further applications can 
Oh! be expected. This leads us to think that it is important and interesting to consider their influences on the physical 
^Jj ■ problems via first reexamining those familiar and standard examples, and then studying more practical systems, 
Ch ' because better precision, higher efficiency, as well as correct physical features are always the aims that physics pursues 
continuously. Here, we focus our attention on a typical example, which shows satisfying results in the calculation 
precision, efficient as well as the corresponding physical features. This implies, in our view, that Wang's scheme of 
perturbation theory is a powerful tool in the perturbative calculation of quantum system. 

In this paper, we intend to study the revisions on transition probability in Ref. which is different from traditional 
one. We find that Wang's revision to the existed expression of transition probability is never trivial, and we illustrate 
this conclusion via calculating this revision for a Hydrogen atom in a constant magnetic field. This example concerns 
Ci I the ground state hyperfine structure of Hydrogen atom whose correction of electron self-energy has been studied (see 
0), so does similar problem about muonium (see [3 and [3). By comparing our result of transition probability with 
traditional one, our above view is verified. After referring to the exact solution of this problem, we will see that the 
improved transition probability given by [l| does show some of its advantages. From this case it is necessary to realize 
that perhaps there are also some other problems that need similar revisions in transition probability. We are sure 
that applications of Wang's scheme of perturbation theory to other problems should not be unimportant and short 
of practical significance, although his scheme just contains the higher order revisions. 

To effectively organize this article, we divide it into the following main parts. Besides Sec. |T]which is an introduction, 
in Sec. |TT]we first introduce the amendment to the transition probability based on the results of [1]. Next, in Sec. IIIII 
we provide the calculation of energy of our example according to the Wang's scheme, then compare it with the exact 
result, and these results are helpful for us later to calculate the transition probability. Then, in Sec. IIVI we will take 
use of Secin] and Sec lIIII to calculate the referred case. Finally, in Sec. |V] we summarize our conclusions and make 
some discussions. 
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II. IMPROVED TRANSITION PROBABILITY 

Let us start with Wang's scheme of perturbation theory denoting the state vector by 
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where |vE'(i)) is the eigenvector of H, and \(P) is the eigenvector of Hq. According to the improved form of the 
perturbed solution of dynamics in [l| , the first order amendment of coefficient of state vector is deduced as 
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4jW = :B^r:^(i-^'"^'*)' (2) 
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where Ej are eigenvalues of Hq with the eigenvector {(f)"'), gj^ (7 7^ (3) are off-diagonal elements of Hi in the 
representation of Hq, that is gj^^ = {(j)^\Hi\(l)p) , and 



1^7/3 = — Ep, (3) 
while E^ is so-called improved form of perturbed energy defined by 

E:,^E., + hJ + + + Gf + • • • . (4) 
Here, hj are diagonal elements of Hi in the representation of Hq, and 



" 0,i:,0,iE0-E,,){E,~E,,){E,~E,,) ^^^{E,-E,,f{Ep-EpS 

where rjpp^ — 1 — Spp^^ rjpp^ = 0, if /3 = /32; 77/3/33 7^ 0, if /? 7^ /32- Thus it leads to so-called improved transition 
probability that is different from the traditional one, that is 



III. REVISION OF PERTURBED ENERGY 

Now we begin to study such a system in which a Hydrogen atom is placed in a uniform constant magnetic field, 
which is in the +z direction in the rectangular coordinate system. Suppose that the atom is at the ground state, thus 
its Hamiltonian (see Ref. [5|] Appendix(Dl)) is 

H = Hq + Hi = B{fieCrez ~ fJ-pCTpz) + WcTe ■ (Tp. (9) 

In the equation above, we change some symbols that are different from the equation (Dl) in Ref. Q. is the 
magnitude of magnetic field. is a coupling constant, and we will give its value later in Sec. IIVI is the electron 
magnetic moment, /ip is the magnetic moment of a proton, and cr — [ux, cFy,az) is the Pauli operator. Since it is well 
known that fj,p <^ fie, we can easily omit the BfipUpz in the H . Therefore, the division of Hamiltonian can be written 
as 

H^Ho+Hi^ Bfie<Jez + W(Te ■ (Tp. (10) 



,7/3 



sin^(w^/3t/2) 
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Here, Hq = WcTe ■ <Tp is the unperturbed part, and Hi = B^eCez is taken as the perturbed part. 

Now we refer to the way proposed by Wang in Ref. [Ij]. First, we should calculate the eigenvalues and eigenvectors 
of Hq. It is easy since we know that Hq is just the coupling of spins of electron and proton. Here we just list the 
results and omit the process of calculating them. 

)i = a{e)a{p) 

Eigenvalue: W =^ Eigenvectors: { (t>2 = + /3(e)a(p)] , (11) 



and 



Here, a and (3 mean: 



mm 



Eigenvalue: — 3W =^ Eigenstate: (j)4 — -^[a{e)l3{p) — f3{e)a{p)] (12) 



V2' 



"=(o) (13) 

It is clear that there is a three-fold degeneracy subspace. 
Considering that 

0-ez|</'l) = |0l), <Jez\(l>2) = \<l>4), CTezlh) = -{h) , <^ez\(l)4) = Ih) (15) 

The H matrix in the representation of Hq can be written as following : 

/W + B^le \ 

W BfiA 

M^~i?Me (1^) 

V B/^e -3W/ 

According to Ref. |l| , we should redivide H into two parts: the diagonal part and the off-diagonal part gi, that 
is, 

^W + B^e \ 

-"o " I W-Bfle 

-3W/ 
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/O 

B/le 



\0 Bfle 



(18) 



It is clear that the degeneracy is completely removed by the redivision skill [l|. 
Again from Ref. Q, we have 



Ep = + Gf + Gf +Gf+--- (19) 



Note that (fl"o)/3/3 ^Ep^- /if. Thus, using the given matrix above, it is easy to obtain that: 

^ = = G^ ^ = 0, (20) 

^(2) _ {Bn^f (3) _ (4) _ {B^lef 

^ -^Vf^' ^2 -U, G2 --j^^ (21) 
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^(2) _ ^(3) _ ^(4) _ „ 



^(2) _ {B^ief ^(3) _ „ ^(4) _ {By.e)'^ 



Thus 



m 
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AW 



{4Wf 



(22) 
(23) 

(24) 
(25) 
(26) 
(27) 



Since this example can be, in fact, solved exactly, the results above obtained via the perturbation theory can be 
compared with the exact solution of energy. In terms of the standard method, we immediately get the exact solution 
of eigenvalues and their corresponding eigenvectors 



Ei 



Ef = W + B^i. 



*1 
1 



El = W - B^ 



Ej ^-W- V4W^' + if^eB)^, =^ *4 = 



^(^^42+^^12)^+4(^^6^)2 

+ *3 = '/>3 
1 



[(^42 + ^^2)02 - 2B/ie(/)4] 



(28) 
(29) 

(30) 



[(CJ42 - '^r2)<^2 - 2B/ie04] (31) 



Here tJ42 = E4 — E2, and 0^X2 = BJ — E^. 

When the magnitude of magnetic field B is very weak, we can expand the exact result of energy, which can be 
compared with the energy worked out by the new way. Assuming that B/ig <C W ^ thus 



while 



El 



~W + 2W{1 



-W + 2W 



1 



{Bf^^^ 

(BMe)' (BfieY 



W 



Eo 



8W^ 8(4W^2)2 

(B/ie)' (Bfie)" 



AW 



[AWf 



Ei = -W -2W [I 



2 \ 2 



= ~W~2W 
« -iW 



1 



4VF2 

{B^Ji,f {B^i,)^ 



m'^ 8(4W^2)2 

{B^Ji,Y , (i?Me)^ 



AW 



{AWf 



Ea 



E. 



(32) 



(33) 



(34) 



Obviously, we see that the perturbed energy calculated by using Wang's scheme [l| is more satisfying because they 
have better precision than the existed method when two kinds of results are compared with the exact results. 
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IV. CALCULATION OF IMPROVED TRANSITION PROBABILITY 

In this section we will consider the transition probability from the state \(f)2) to state |04). First, we should calculate 
its exact result. Obviously the transition probability is: 



P^(2^4) = |((^4|e 



-iHt\ 



(35) 



To work out this equation, we take use of the completeness of the Hilbert space {|^'i), |^'2), l^a), 1^4)} given above, 
which means: 



P^(2 ^ 4) 



71.72 



= 1 \'^4||^7iM^7i 



!'4||*2)(*2||02)e-''^"* + (04||*4)(«'4||02)e-*^"* 



■ 2 sin^(ijj|2i/2) 



(36) 



K2/2)^ 



where — EJ — = —2^4:W'^ + {fieBY- We also write down the result given by the traditional perturbation 
theory: 



P(2 ^ 4) - {^i,B) 



, sin^(tJ42t/2) 

(^42/2)2 



where W42 = E4 — E2 = ~AW. and the result given by Eq.® is: 



(37) 



(38) 



/r) \2 / \4 

where W42 = E4 — E2 = —AW j^^^^ h 2^777^. Since Eqs. ([36]), (jST]) and ((38)) have some common factors such as 



2W 



(^eB)^ and o;?,! order to effectively compare them we can let the three equations plus a common factor: . 

[^J'eBy 

Besides, to calculate we should now add h at the right place (since we have let h = I previously). That is to say, we 
just need to compare , Pj, and P, which are: 



= P^(2 ^ 4) X 



(cJ42/2)2 _ sm'{uji2t/2h) 



(39) 



4VF2 



P7 = Pj(2^4)x 



(c^42/2)^ 

{f^eBr 



sm'^{LU42t/2h) = sin^ 



2w + i^-iM-]t/n 



{my 



(40) 



P = P{2~^ 4) 



(^42/2)^ 

(/^eP)2 



: sm 



'{LU42t/2h) = sin^ i2Wt/h) 



(41) 



Next, we will replace the parameters such as W and B with concrete data. From the latest data book we get 
the value of = 9.28476412 x 10"2''J-T'^ Now we explain the W in Eq. dTO]). Comparing with the equation (Dl) 
of Ref. [5], Eq. (fTUl) here is written in the form of Pauli operator but not of spin operator which is used in (Dl) of 
0. Therefore we can easily know that: 



W 



2tt 



(42) 



where the Ai'h is the so-called ground-state hyperfine frequency, and h is the Plank constant. By referring to Q we 
get the experimental value of A^h = 1420.4057517667 MHz. Plus the h = 6.6260693 x 10"34gY.g (ggg xABEL XXV 
of 01), we finally work out the value of W: 



W = 1.46858145124* lO^'^eV 



(43) 
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A. COMPARISON OF TRANSITION PROBABILITY AS TIME GOES 

Now we want to compare Eqs. ((39|) . (|4ip and (|40)) along with the change of time t and magnetic field B respectively. 
First, we compare them when time is changing, and we set the magnitude of magnetic field as B — which can 

be easily realized in laboratory. With all the data we get, we now rewrite the expression of P'^, Pi and P: 

_ sin^ (4.46320474159 x 10^ x t) 
1 + 3.89282563044 x IQ-^ 
Pi sin2(4.46320474158x 10^ X t) (44) 

P ^ sin2(4.46233627125 x lO^ x t) 

To see more clearly, we use mathematica tool to show their differences (here we retain 12 effective figures). Note 
that the function sin(t) is an oscillating one, and that the period of is so small that the three functions of 
will oscillate intensively if the parameter t varies in a large range. This make us have to compare them within a small 
time span. The graphs bellow show four small-time-spans, in which we can see clearly their nuances(see Fig.l). 



Pic-A Pic-B 




FIG. 1: The four graphs represent four small-time-spans, red curve represents the P"^, green the Pi, and blue the P. We do 
not mark the ticks because the time span is so small that the Mathematica tool can not even mark them clearly,yet this does 
not imped us to see the differences among the three curves of the four graphs above. The Pic-A shows while t is around 10~^s, 
P^ , Pi and P coincide with each other. Pic-B shows while t is around Is, P^ and Pi still coincide with each other, but P 
begins to deviate. Pic-C shows when t is around 6s, Pi begins to stray away from P^ , while P also deviates from the former 
two. Pic-D shows when t is around 27.7s, obvious deviations between P^ and Pi occur, and P still shows irregular deviation. 

FIG.l tells us that when time t is very small, say 10~^s, the transition probability of traditional way and new 
way both coincide with the exact solution. However, when time evolves to a relatively large value, traditional result 
show irregular deviation from the exact one, while new result still matches well with it. Of course, new result also 
will deviate from the exact one when time t is large enough. Therefore we can say that, when the magnitude of 
magnetic field is given, traditional result is useful only when time t is small enough, but new one will be exact in a 
large time-span. As a result, we clearly see the advantage of new way of perturbation theory. Next, we will compare 
the three curves when time "t" is set and magnetic field is changing. 

B. COMPARISON OF TRANSITION PROBABILITY AS MAGNETIC FIELD CHANGES 

At the given time we choose the time t = Is. We also choose four small magnetic field spans, and FIG. 2 represents 
the comparison of three curves: P^, P/ and P. Here we write down the expression of them: 
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J, _ sin2(Vl. 991244499780254521326 x IQis + 7.751567612132237407735 x lO^i x B^) 

1 + 389.2825623868 x 

Pi = sin2(4.462336271259 x 10^ + 8. 685548489539 X 10" X S2 - 8.452831429358 X 10" X (45) 
P = sin2(4.462336271259 x 10^) 



Pic-1 Pic-2 




FIG. 2: In the graph above, the red curve represents , green the Pj and blue the P. Curve P is a straight line which is 
independent to B. Pic-1 shows when B is around 10~*T, P^ and Pj coincide with each other. Pic-2 shows when B is around 
1.29 X 10"^T, nuance between P^ and Pj appears. Pic-3 shows when B is around 1.21 x 10 ^T, there is obvious distinction 
between them. Pic-4 shows when B is around 0.036T, huge difference occurs, and even their oscillation amplitudes are different. 

In this part, we again see the differences between traditional result and new one. In fact, in this case the phase 
factor of traditional transition probability has nothing to do with magnetic field, which results in the straight line 
in FIG. 2. However, as we can see from the exact solution, the transition probability oscillate intensively with the 
change of B, and this reveals the defect of traditional way. Yet, new way in [T] provide us with a relatively satisfying 
expression, which also vibrates with magnetic field "-B" . This new expression matches well with exact solution when 
"B" varies in a certain range. Yet, we should notice that there is also some flaw in the new result, which can be seen 
from the Pic-4 in FIG. 2. Such flaw results from the difference of denominators between Eqs. ([5]) and (|36p . which will 
influence their vibration amplitudes. 



V. DISCUSSION AND CONCLUSION 

Now we will come to a conclusion. In Ref. [l[, the author gives a improved form of perturbed solution of dynamics, it 
is interesting to study its influences on some physics results. As we can see, the phase factors correlating to dynamical 
behavior in the expression of state vector are changed, and the variations absorb the partial contributions from higher 
order approximations. This leads us to consider the change of transition probability. Since abstract discussion of the 
transition probability is both hard and has no common meaning, we should find concrete cases to discuss. 

Main work we do in this paper is to calculate the improved transition probability for a Hydrogen atom in a constant 
magnetic field. In fact, this example we choose is appropriate, because we need to choose one which can be given the 
exact solution so that we can compare it with the result of traditional and Wang's scheme of perturbation theory. The 
comparison we give in Sec. [IV] successfully indicate that Wang's scheme does show its advantages beyond traditional 
one. Perhaps someone would argue that this example is too simple to be meaningful enough, however, we should 
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note that such simple case do suggest that some expressions such as transition probabihty in other cases might need 
amended, or at least need perfected. Maybe such amendment does not have huge influence in some systems, but to 
extend this assertion to all problems is both too early and too cursory. 

In fact, Ref. [ij uses several skillful ways to include the high-order approximation in the expression of low-order 
one, and this procedure results in a tidy expression of perturbed solution of dynamics. It is a physical reason why 
Wang's scheme can improve the existed some conclusions. Finally, since we also see the flaw shown in Pic-4 of Fig. 2, 
perhaps there will still be further amendment to it in the future. 
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